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Preface 


The  purpose  of  this  study  is  to  analyze  the  predictions  made  by  the  Space  Surveillance  Center 
(SSC)  concerning  the  impact  time  and  location  of  decaying  satellites.  Because  of  my  prior 
experience  m  the  SSC  and  because  I  believe  there  to  be  a  general  tack  of  confidence  by  the  public  in 
these  predictions,  I  was  very  interested  in  determining  how  accurate  these  predictions  really  are. 

There  are  two  specific  goals  for  this  thesis.  The  first  goal  is  to  directly  compare  the  accuracy  of 
the  Tracking  and  Impact  Predictions  (TIPs)  made  by  the  SSC  during  the  four  years  from  1987  to 
1990  to  the  accuracy  level  asserted  by  the  SSC.  The  second  goal  is  to  determme  if,  based  on  these 
fmdings,  it  would  be  a<h’antageous  for  the  SSC  to  prepare  its  OPR£P*3  report  earUer  than  is 
currently  being  done.  The  OPR£P«3  report  is  used  to  notify  the  proper  authorities  of  possible 
impact  by  a  deca)ing  object  near  or  witlun  the  Soviet  border. 

In  my  attempt  to  perform  the  analysis,  I  received  a  great  deal  of  help  from  seveial  people  to 
whom  I  would  like  to  express  my  sincere  gratitude  I  would  first  like  to  state  that  I  am  deeply 
indebted  to  my  thesis  advisor.  Professor  T.  S.  Kelso,  for  his  constant  patience,  insight,  and  direction. 
Second,  I  would  like  to  thank  Dr.  William  E.  Wiescl  for  sharing  his  knowledge  of  the  programs  used 
by  the  SSC.  Third,  I  would  like  to  thank  Professor  Daniel  E.  Reynolds  for  his  undying  enthusiasm 
and  encouragement  as  well  as  his  statistical  ei^rtise.  He  has  that  precious  ability  to  make  what  can 
sometimes  seem  very  <Ufficult  and  time  consuming  actually  fun.  Finally,  I  would  like  to  thank  my 
fianc6  Charles  J.  Martin,  Jr.,  the  most  important  person  in  my  life,  for  his  constant  support  and 
understanding  during  those  many  months  when  I  was  glued  to  either  the  microfiche  reader  or  the 
computer  saeen. 


Susanne  V.  Lefebvre 
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Abstract 

This  study  analyzed  the  accuracy  of  the  early  Tracking  and  Impact  Predictions  (TIPs)  made  by  the 
United  States  Space  Command's  Space  Surveillance  Center  (SSC)  during  the  1987>1990  time  period. 
The  fmal  prediction  data  was  first  compared  to  visual  observation  data  (Vis  Obs)  for  those  objects  for 
which  Vis  Ob  data  was  available.  The  early  TIP  decay  predictions  which  include  the  7-day  through  3- 
hour  predictions  were  then  compared  to  the  final  prediction  for  each  TIP  object.  For  each  TIP  run,  the 
time  error  was  then  calculated  as  a  percentage  of  the  time  until  decay.  The  results  were  then  compared 
to  the  accuracy  level  asserted  by  the  SSC 

The  results  of  this  study  indicate  that  the  accuracy  of  the  decay  predictions  is  usually,  but  not  always, 
withm  the  ±20  percent  accuracy  level  asserted  by  the  SSC.  The  results  also  suggest  the  existence  of  a 
positive  bias  indicating  that  the  early  TIP  decay  predictions  are  routinely  late  relatKY  to  the  final  decay 
prediction. 

An  attempt  was  then  made  to  model  cut  some  of  the  positive  bias  found  is  the  TIP  decay  prediction 
data  using  multiple  linear  regression.  Six  regression  models  were  found  which,  if  incorporated  into  the 
current  SSC  TIP  decay  procedures,  would  allow  the  SSC  to  predict  the  final  decay  prediction  time  with 
substantially  less  error. 

’Hiis  study  also  analyzed  ti  s  current  SSC  OPREP-3  report  initiation  procedure.  It  was  determined 
that  the  dedsion  to  initiate  the  report  should  remain  at  the  6-hour  point  but  that  one  of  the  regression 
models  mentioned  above,  namely  Model  6,  should  be  used  in  conjunction  with  the  TIP  decay  prediction 
data  to  obtain  a  better  estimate  of  (he  final  decay  time.  It  was  also  determined  that  the  ±15  minute 
error  window  currently  used  as  a  guide  for  determining  the  necessity  of  the  OPREP-3  report  should  be 
widened  to  approximately  ±40  minutes  to  more  accurately  account  for  error  extremes  in  the  TIP  decay 
prediction  calculations. 


AN  assessment  OF  TRACKING  AND 


IMPACT  PREDICTIONS 


I.  Introduction 

Backffvund 

The  mam  mia&ioo  of  the  Space  SurveiUance  Center  (SSC),  located  at  Cheyenne  Mountain  Air 
Force  Base,  Colorado,  is  to  detect,  track,  identify,  and  catalog  all  man*madc  space  objects  in  support 
of  the  United  States  Space  Command  (USSPACECOM).  One  of  the  responsibilities  of  the  SSC  is 
to  provide  Tracking  and  Impact  Prediction  (TIP)  for  those  objects  that  are  within  fifteen  days  of 
decay  and  which  have  a  high  probability  of  impaamg  the  earth.  The  purpose  of  TIP  processing  is  to 
ensure  ‘that  the  Missile  Warning  Center  and  the  USSR  (when  applicable)  will  not  mistake  a  TIP 
object  for  an  incoming  RV  (re*eotry  vehicle)’  (6  Sec  9,13)  and  ”to  determine  which  country  is  hable 
for  damages  resulting  from  a  satellite’s  impact’  (6  Sec  9,13).  TIP  processing  uses  special 
perturbations  theory  to  predict  the  decay  time  and  location  of  all  pa^oads,  rocket  bodies,  platforms, 
and  debris  with  a  greater  than  five  percent  chanw  of  surviving  re*entry. 

The  impaa  predictions  are  updated  at  speaflc  time  intervals  throughout  the  TIP  objea’s  decay. 
Each  new  prediction,  called  a  TIP  mn,  includes  a  decay  time,  a  decay  location,  and  a  ground  trace. 
The  ground  trace  is  a  map  upon  which  the  decaying  satellite’s  flight  path  is  drawn  and  is  used  to 
determine  \riiether  the  USSR  should  be  notified.  If  deemed  necessary,  the  information  will  be 
passed  up>channel  to  the  proper  authonties  in  an  OPREP'3  report.  The  more  accurate  the  early 
predictions,  the  sooner  the  proper  authorities  can  be  notified.  The  author,  having  worked  as  an 
orbital  analyst  in  the  SSC,  has  noted  a  genera!  lade  of  confidence  in  the  accuracy  of  the  earlier 
predictions  and  the  use,  therefore,  of  the  later  predictions  for  such  decisions  as  the  necessity  of 
OPREP>3  report  initiation. 
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Reseanh  Objective 

The  purpose  of  this  research  project  is  two-fold:  (1)  to  compare  the  accuracy  of  the  TIP  decay 
predictions  for  the  1987-1990  time  period  to  the  accura^  level  asserted  by  the  Space  Surveillance 
Center  and  (2)  to  determine  if  it  would  be  advantageous  for  the  SSC  to  initiate  the  OPREP-3  report 
earlier  than  is  currently  done. 

Overview 

The  remainder  of  this  study  includes  a  literature  review  of  orbital  mechanics>  perturbative  forces, 
differential  corrections,  SSC  satellite  decay  processmg,  and  multiple  linear  regression  (Chapter  II), 
data  collection  and  description  (Chapter  III),  methodology  used  to  conduct  the  study  (Chapter  IV), 
results  of  the  analysis  (Chapter  V),  and  finally,  conclusions  and  recommendations  (Chapter  VI). 
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II.  Review  of  Liunture  and  Background  Development 


Introduction 

rfac  following  pages  will  review  literature  oa  five  topics  pertinent  to  this  research  proposal.  The 
intent  of  this  section  is  to  expand  the  background  information  presented  earher.  The  speofir  topics 
discussed  are  orbital  mechanics,  perturbative  forces,  differential  corrections,  satelhte  decay 
proces^ng  as  performed  by  the  Space  Surveillance  Center  (SSC),  and  multiple  linear  regression. 

Orbital  Mechanics 

This  section  discusses  the  concepts  and  definitions  of  orbital  mechanics  necessary  for  a  basic 
understanding  of  the  orbital  motion  of  an  artiffcial  satellite  about  the  earth.  Five  independent 
parameters,  called  orbual  elements^  are  required  to  completely  describe  the  size,  shape,  and 
orientation  of  an  orbit.  To  further  pinpoint  the  position  of  a  satellite  along  its  orbit  at  a  particular 
time,  a  sixth  parameter  is  required.  Together,  these  six  parameters  form  an  orbital  element  set  and 
allow  one  satellite  to  be  distinguished  from  thousands  of  other  satellites  in  earth  orbit.  One  type  of 
orbital  element  set  is  the  Keplerian  element  set.  it  includes  the  semi^major  axis,  eccentricity, 
inclination,  r^t  ascension  of  the  ascending  node,  argument  of  perigee,  true  anomaly,  and  epoch 
time  (3:58, 6:Sec  2,5)  (see  Figure  1).  The  semi^inajor  axis  (a)  is  one-half  the  longest  diameter  of  a 
satellite’s  orbit  and  is  used  to  describe  the  rize  of  the  orbit.  The  longer  the  semi-major  axis,  the 
larger  the  orbit.  The  semi-major  axis  can  also  be  used  to  determine  the  time  required  for  the 
satellite  to  complete  one  revolution.  This  quantity  is  known  as  the  satellite’s  orbital  period.  For 
convenience,  the  SSC  has  classified  artificial  earth  satellites  into  three  general  categories:  deep- 
space  satellites,  near-earth  satellites,  and  decaying  satellites.  Decp-space  satellites  are  those 
satellites  with  orbital  periods  greater  than  or  equal  to  225  minutes.  Near-earth  satellites  are  those 
satellites  with  orbital  periods  between  87.5  minutes  and  225  minutes.  Decaying  satellites  are  those 
satellites  vrith  orbital  periods  equal  to  or  less  than  875  minutes  (10). 
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Figure  1.  Orbital  Elements  (Reprinted  from  7iec  2^) 

Eccentricity  (c)  is  a  parameter  used  to  describe  the  shape  of  a  satellite’s  orbit.  Eanh  sateilites 
can  have  only  circular  or  elliptical  orbits.  The  eccentricity  of  a  perfectly  circiilar  orbit  is  equal  to 
zero,  while  that  of  an  elliptical  orbit  is  between  zero  and  one.  The  closer  the  eccentricity  is  to  one, 
the  more  elongated  the  orbit. 

Inclination  (i)  is  one  of  two  parameters  used  to  locate  the  orbital  plane  in  space.  It  is  the 
angular  measurement  between  the  equatorial  plane  and  the  orbital  plane  and  is  measured  in  a 
counterdocicwise  direction  at  the  ascending  node.  Simply  put,  the  inefination  is  the  angle  between 
the  orbit’s  angular  momentum  vector  and  the  earth’s.  The  indinalion  determines  whether  an  orbit 
is  prograde  or  retrograde,  polar  or  equatorial  If  the  inclination  is  greater  than  or  equal  to  0*  but 
less  than  9C*,  the  orbit  is  prograde.  If  the  inclination  is  greater  than  90*  but  less  rl*a"  or  equal  to 
180*,  the  orbit  is  retrograde.  If  the  inclination  is  near  90*,  the  orbit  is  polar.  If  the  indinalion  is 
equal  to  0*,  the  orbit  is  equatorial  (6:Sec  2^). 
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Right  ascension  of  the  ascending  node  (a)  is  the  second  parameter  used  to  locate  the  orbital 
plane  in  space.  It  is  an  angular  measurement  in  the  equatorial  plane  from  the  vernal  equinox  to  the 
ascending  node.  The  vernal  equinox  is  an  imaginary  line  drawn  in  inertial  space  from  the  center  of 
the  earth,  through  the  equator,  towards  the  sun  at  the  bepnning  of  spring.  This  direction  is 
routinely  referred  to  as  the  first  point  of  Anes  (T)  (see  lugure  2).  The  ascending  node  is  the  point 
at  which  a  satellite's  giound  trace  intersects  the  equator  as  the  satellite  travels  from  the  southern 
hemisphere  into  the  northern  hemi'.-'  •'e*  (6:Sec  2,6). 


TO  THE  FIRST 
POINT  OP  ARIES  (T) 


Figure  2.  The  Vernal  Equinox  (Reprinted  from  7;Sec  2,28) 

Argument  of  perigee  (u)  is  an  angle  used  to  orient  the  orbit  within  the  orbital  plane.  It  b  the 
angle  that  b  swept  out  by  the  satellite  as  it  traveb  from  the  ascending  node  to  Ih:  perigee  point 
Perigee  b  the  position  in  an  earth  satellite’s  orbit  that  b  closest  to  the  center  of  the  earth.  It  should 
be  noted  that  perigee  b  not  deSned  for  a  perfectly  circular  orbit  since  all  points  in  a  drcular  orbit 
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are  an  equal  distance  from  the  center  of  the  earth.  It  should  also  be  noted  that  in  practice,  there  is 
no  such  thing  as  a  perfectly  circular  orbit. 

True  Anomaly  (v)  is  the  angular  measurement  in  the  direction  of  the  satellite's  motion  along  its 
orbital  plane  from  perigee  to  the  satellite’s  position  at  epoch  time”  (6:Sect  2,8).  Epoch  time  can  be 
any  arbitrary  moment  m  time  and  is  used  as  a  reference  pomt.  In  the  Space  Surveillance  Center, 
the  epoch  time  is  usually  set  to  the  time  the  ol]^ect  was  last  observed  (1).  True  anomaly  and  epoch 
time  go  together.  One  is  of  no  use  without  the  <»her.  The  SSC  does  not  directly  use  true  anomaly, 
however.  Instead,  it  uses  a  parameter  known  as  neon  onomn/y,  (M).  To  derive  mean  anomaly, 
another  angle,  eccentric  anamafy^  (E)  is  first  determined  usin^  the  geometric  relationship  shown  in 
Figure  3  The  eccentric  anomaly  is  then  used  to  derive  the  mean  anomaly  uriog  Equation  1  (3:183* 
184,  6:23*9). 

Af=£-esm(£)  U) 

where 

E  is  the  eccentric  anomaly 
e  is  the  eccentricity 

It  should  be  noted  that  unlike  true  and  eccentric  anomaly,  mean  anomaly  is  not  an  angle  but  a 
mathematical  relationship.  The  use  of  mean  anomaly  instead  of  true  anomaly  allows  a  satellite’s 
orbital  path  to  be  modeled  on  a  circle  instead  of  an  ellipse.  In  this  way,  the  satellite  is  modeled  as 
moving  at  a  constant  velocity.  If  an  elliptical  model  were  used,  a  more  complicated  mathematical 
model  would  be  needed  to  account  for  the  gradual  increase  in  velocity  of  a  satellite  as  it  approaches 
perigee  and  the  gradual  decrease  in  velocity  as  it  approaches  apogee.  Mean  anomaly  is  used, 
therefore,  to  simplify  the  problem  of  predicting  where  a  satellite  will  be  in  the  future  (6:23). 

Perturbative  Forces 

There  are  three  major  perturbative  forces  that  affect  the  orlrit  of  most  artificial  earth  satellites. 
The  magnitude  of  these  forces  depends  on  the  satellite’s  size,  shape,  mass,  and  orbit.  The  three 
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Figure  3.  Eccentric  Anomaly  (Reprinted  from  3:183) 


forces  arc  the  direct  result  of  atmospheric  drag,  the  nonsphcriciJ  shape  of  the  earth,  and  third-body 
gravitational  attraction  Together,  these  perturbattons  not  only  change  a  satellite’s  orbit,  but  ran  also 
cause  the  satellite  to  plummet  to  earth  or  decay  (16:165). 

Atmospheric  drag  is  caused  by  the  Metion  between  the  molecules  in  the  earth’s  upper 
atmosphere  and  the  surface  of  the  satellite.  Atmospheric  drag  acts  m  a  dircaion  opposite  to  the 
velocity  of  the  satellite  relative  to  the  atmosphere.  Because  the  near-earth  environment  is 
characterized  by  a  higher  atmospheric  density  than  the  deep-space  environment,  near-earth  satellites 
experience  greater  atmospheric  drag  than  do  deep-space  satellites.  When  referring  to  the 
atmospheric  drag  on  a  satellite,  a  dimensionless  quantity  {C^),  called  the  drag  coefficient,  is  used. 
The  drag  coefficient  is  related  to  the  shape  of  the  satellite  (3:423-424). 

The  remaining  two  perturbative  forces,  those  resulting  from  the  nonspherical  shape  of  the  earth 
and  third-body  attractions,  are  less  sigmOcont  for  near-earth  satellites  than  for  deep-space  satellites 
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The  nonspherical  shape  of  the  earth  is  the  result  of  the  fact  that  the  earth  has  a  noouniform  shape 
and  density.  This  phenomenon  results  in  an  asymmetrical  mass  distribution  that  causes  the  earth’s 
gravitational  force  on  the  satellite  to  be  directed  slightly  away  from  the  center  of  the  earth.  The 
third-body  gravitational  effects  are  due  to  the  fact  that  the  earth  is  not  the  only  gravitational  force 
acting  on  a  satellite’s  orbit.  Other  planetary  bodies  also  exert  a  gravitational  pull  on  the  sateUite. 
These  third-body  effects  are  stronger  for  closer  bodies  such  as  the  sun  and  moon  (6:Sec  64)* 

Differential  Comctions 

The  Space  Surveillance  Center  (SSC)  uses  differential  corrections  to  update  its  element  sets. 
Element  sets,  or  elsets^  change  daily,  mainly  due  to  perturbations  such  as  those  previously  discussed. 
These  d^y  changes  must  be  incorporated  into  the  element  sets  and  are  done  so  through  the  use  of 
differential  corrections.  Satellite  observations  sent  by  the  Space  Sensor  Network  (SSN)  and  received 
by  the  SSC  do  not  exactly  fit  the  position  predicted  by  the  element  sets.  These  deviations  are  caused 
by  changes  in  a  satellite’s  orbital  elements  due  to  perturbations  which  were  not  accounted  for  by  the 
model,  by  the  fact  that  the  element  set  and  the  associated  perturbation  modeling  cannot  predict  the 
exact  position  of  the  satellite,  and  by  *each  sensor’s  error  (bias)  in  measuring  the  exact  position  of  a 
satellite"  (6:6,6).  Differential  corrections  are  used  to  update  element  sets  by  attempting  to  eliminate 
any  reriduals  caused  by  changes  to  orbital  elements.  To  minimize  errors,  a  least  squares  fit  is  used. 
Simply  put,  "differential  corrections  attempt  to  mathematically  fit  the  best  ellipse  (with  associated 
perturbations)  to  the  observations  of  a  particular  satellite"  (6:6,6). 

The  two  main  perturbation  models  used  by  the  SSC  when  performing  its  differential  conectioos 
are  the  General  Pertutbadons  (GP)  model  and  the  Special  Perturbadons  (SP)  model  For  both  GP 
and  SP  elsets,  an  elliptical  orbit  is  found  such  that  when  the  respective  perturbations  are  added  on, 
the  result  will  represent  the  true  path  of  the  satellite  as  closely  as  possible.  Because  SP  modeling 
accounts  for  more  perturbations  than  does  GP  modeling,  the  final  path  described  by  an  SP  elset  will 
fit  the  satellite’s  true  path  better  than  that  described  by  a  GP  elset. 
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General  Perturbations  Modeling.  The  GP  model  uses  analytical  equations  to  incorporate  latitude- 
dependent  perturbations.  These  perturbations  include  the  imbalanced  mass  distribution  of  the  earth 
due  to  its  nonspberical  shape,  the  greater  mass  present  in  the  southern  hemisphere  of  the  earth,  and 
other  observed  mass  anomalies.  Atmospheric  effects  zxt  accounted  for  by  using  the  BSTAR  drag 
term  (6.Sec  6,9).  BSTAR  is  derived  from  the  atmospheric  model  and  is  a  ballistic  coefficient  for 
atmospheric  drag  (20  5).  The  units  of  BSTAR  are  1/ER,  where  ER  is  the  earth’s  average  equatorial 
radius  (19:123). 

Special  Pertu/bations  Modeling  The  SP  model  emplt^  numerical  methods,  rather  than 
analytical  methods,  to  incorporate  latitude-dependent  perturbations  (zonal),  longitude-dependent 
perturbations  (sectoral),  and  both  latitude  and  longitude-dependent  perturbations  (tesseral). 
Gravitational  effects  are  accounted  for  by  breaking  the  surface  of  the  earth  down  into  small  grids  to 
allow  for  better  resolution.  Gravitational  effects  of  the  sun,  moon,  and  planets  can  also  be  modeled. 

Occasionally,  satellites  will  pass  throu^  the  same  area  and  experience  the  same  perturbative 
forces.  When  this  occurs,  the  satellites  encounter  resonance  effects.  The  SP  model  is  also  capable 
of  accounting  for  this  phenomenon. 

To  model  atmospheric  effects,  SP  modeUi^  uses  the  Jacchia-Nicolet  model,  which  is  more 
sophisticated  than  the  atmospheric  model  used  for  GP  updates.  The  Jacchia-Nicolet  model  accounts 
for  the  earth’s  diurnal  bulge,  solar  activity,  geomagnetic  activity,  and  semiannual  variation.  The 
diurnal  bulge  is  due  to  the  fact  that  the  sunlit  side  of  the  earth  is  warmer  than  the  dark  side,  causing 
the  atmosphere  on  the  sunht  ride  to  expand  into  space.  This  expansion  results  in  a  varying 
atmospheric  density  based  on  altitude  and  sun  angle.  The  level  of  solar  activity  is  also  incorporated 
into  the  model  to  account  for  the  inaease  in  the  number  of  charged  particles  ejected  from  the  sun 
at  higher  levels  of  solar  activity.  Geomagnetic  activity  is  accounted  for  1^  incorporating  the  earth’s 
magnetic  field  strength,  measured  in  Ap  levels,  wbidi  is  also  affected  by  solar  activity.  Semiannual 
variation  refers  to  measured  harmonic  variations  in  the  atmosphere  due  to  seasonal  changes  (6:Sec 
6,9-10).  It  should  also  be  noted,  however,  that  there  are  limits  to  the  predictive  ability  of  the 
atmospheric  model.  One  major  source  of  TIP  prediction  error,  for  example,  is  related  to  our 


9 


current  mability  to  accurately  predict  future  solar  activity  and  its  subsequent  ejects  on  our 
atmosphere.  Currently,  the  best  we  can  do  is  predict  changes  in  the  atmosphere  which  would  occur 
about  two  days  after  an  observed  solar  event.  Because  the  atmospheric  model  cannot  accurately 
account  for  hiture  solar  activity,  errors  in  the  predicted  decay  tunes  can  be  substantial. 

Satellite  Decay  Processing 

This  section  discusses  the  methods  used  by  the  Space  Surveillance  Center  to  predict  the  time 
and  location  of  decaying  satellites  The  SSC  processes  two  types  of  satellite  decays.  The  first  type, 
Normal  Decays^  refers  to  ail  objects  except  pa^oads,  rocket  bodies,  platforms,  and  significant  debris 
(those  with  radar  aoss  sections  greater  than  one  square  meter)  which  have  a  predicted  decay  date 
within  Uurty  days.  These  objects  have  less  th^  a  five  percent  chance  of  surviving  re-entry.  The 
second  type  of  satellite  decay  is  known  as  Tracking  and  Impact  Prediction  (TIP).  TIP  decays  include 
all  pa^oads,  rocket  bodies,  platforms,  and  si^iificant  debris  vriuch  have  a  predicted  decay  date 
within  fifteen  days.  These  objects  have  a  greater  than  five  percent  chance  of  surviving  re>entry  and, 
therefore,  a  greater  likelihood  of  causing  damage  upon  impact.  For  both  types  of  decays,  the  SSC 
will  manually  update  Ue  orbital  elements  of  the  decaying  objects  before  a  decay  prediction  is  made. 
The  processing  of  these  two  types  of  decays  ihlfers  in  the  sophistication  of  the  perturbative  model 
used  for  the  clement  set  updates,  the  frequency  of  updates,  and  the  programs  used  to  compute  the 
predicted  decay  parameters. 

Normal  Decay  Processing.  The  SSC  runs  a  program  called  DECAYX  to  identify  new  normal 
decays  and  compute  the  predicted  decay  dates  of  all  normal  decays.  The  program  is  run  daily.  The 
clement  sets  of  these  o(^ects  are  then  updated  using  general  perturbations  techniques.  DECAYX  is 
rerun  once  all  the  element  sets  of  the  normal  decays  have  been  updated.  The  new  decay  dates  are 
then  logged  for  future  reference.  If  a  normal  decay  is  no  longer  tracked  the  Space  Sensor 
Network,  and  it  is  past  its  most  current  predicted  decay  date,  it  is  assumed  to  have  burned  in  and  is 
decoyed  from  the  SSC’s  active  satellite  database. 
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TIP  Decay  Processing.  Because  TIP  decays  have  a  higher  probability  of  surviving  re-entry,  their 
element  sets  are  updated  using  special  perturbations  techniques.  TIP  procesUng  is  more  rigorous 
than  normal  decay  processing  due  to  the  higher  probability  of  a  TIP  object  surviving  re-entry.  Once 
a  satellite  has  been  identified  as  a  TIP  object,  SP  element  set  updates  and  decay  predictions  are 
performed  at  specific  points  in  time:  seven-to-ten  days  before  re-entry,  four  days  before  re-entry, 
two  days  before  re-entry,  one  day  before  re-entry,  twelve  hours  before  re-entry,  six  hours  before  re¬ 
entry,  two-to-three  hours  before  re-entry,  and  just  after  re-entry.  The  programs  used  by  the  SSC  to 
update  the  orbital  clement  sets  and  then  compute  the  decay  parameters  for  TIP  objects  are  called 
MANDC  and  TIPX,  respectively.  Each  updated  decay  prediction  includes  a  decay  latitude  and 
longitude,  a  decay  time,  and  a  decay  time  error  window.  A  ground  trace  is  also  produced  for  all  but 
the  seven-to-ten  day  run. 

The  difference  between  one  TIP  run  and  the  next  successive  TIP  run  essentially  lies  in  the 
obspan  used  for  the  SP  MANDC  update.  The  obspan  is  the  time  period  from  which  the 
observations  (which  are  used  to  update  the  SP  element  set)  are  chosen.  The  obspans  used  for  the 
vanous  TIP  runs  are  shown  below  in  Table  1.  The  seven-to-ten  day  TIP  run,  for  example,  uses  a 
two-day  obspan  which  means  that  the  SP  MANDC  will  use  the  last  two  days  of  observations  in 


TABLE  I 

TIP  RUN  OBSPANS  (5Atch  1,2) 


1  TIP  RUN 

OBSPAN  (days) 

7-10  day 

2.00 

4  day 

1.40 

2  day 

125 

1  day 

1.00 

12  hr 

OBO 

6  hr 

0.67  1 

2-3  hr 

050  1 

Final 

050  1 
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Its  update.  The  number  of  observations  within  a  particular  obspan  will  vary  from  one  object  to 
another. 


As  a  TIP  object  approaches  its  impact  pointy  the  Space  Sensor  Network  is  periodically  instructed 
to  inaease  the  tasking  for  the  object.  The  tasking  is  a  means  by  which  observational  requirements 
are  set  for  each  satellite.  The  tasking  procedure  is  necessary  because  there  are  only  a  few  sensors 
with  which  to  track  the  large  number  of  orbiting  satellites.  Tasking  ensures  that  the  proper  number 
and  dispersion  of  observations  will  be  obtained  on  each  satellite  and  that  observations  on  high* 
interest  satellites  will  be  obtmned  and  forwarded  to  the  SSC  on  a  priority  basis.  Tasking  codes 
include  a  category  and  a  suffix.  The  category  sets  the  priority  (or  importance)  of  obtaining  and 
transmittir^  the  observations.  The  categories  range  from  CAT  1  to  CAT  3  with  CAT  1  being  used 
for  events  of  highest  priority  and  CAT  3  for  routine  near-earth  satellites  (6:Sec  3,11)-  The  suffix 
delines  the  amount  of  observational  data  required  for  the  satellites.  The  suffixes  used  and  their 
mcanii^  are  listed  in  Table  2.  The  tasking  requirements  for  TIP  decay  processing  is  shown  in 
Tabic  3. 


TABLE  2 

TASKING  SUFFIXES  (6:Scc  3,11-12) 


SUFFIX 

AMOUNT  OF  OBSERVATIONAL  DATA 

A 

Maximum  data  on  all  avmlabie  passes. 

B 

1-1/2  minutes  of  data  on  aU  available  passes. 

D 

Three  observations  evenly  spaced  throughout  the  track  on  all  available 
passes.  (Phased  Arra>'s) 

H 

Three  observations  evenfy  spaced  Ihrougbout  the  track  on  all  available 
passes.  (Mechanical  Trackeis) 

M 

Used  for  deep-space  maneuverable  satellites.  Up  to  20  observations 
are  required  for  a  suspeaed  or  detected  orbit  deviation. 

S 

As  specified  in  a  message  by  the  SSC  for  n»r-eatth  satellites. 

Requires  an  in-track  or  cross-track  search  for  deep-space  satellites. 

T 

Maximum  data  on  all  avjulable  passes  (used  for  later  stages  of  TIP 
decay). 
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Note  that  each  successive  TIP  run  uses  a  smaller  obspan.  Although  the  obspao  is  smaller,  the 


number  of  observations  in  that  obspan  will  not  necessarily  be  less  due  to  the  change  m  tasking 
requirements.  The  observations  will  also  become  closer  in  time  to  the  actual  decay  time 
and  will,  therefore,  more  accurately  reflect  the  current  decay  path  of  the  object. 

It  should  also  be  noted  that  intermediate  updates  are  often  run  on  most  of  the  TIP  objects.  The 
extra  updates  are  mostly  nm  between  the  7*day  run  and  the  l-day  run,  but  can 


TABLE  3 

TIP  TASKING  (5Alch  1,1) 


DAYS  FROM  DECAY 

TASKING 

>  2 

2B 

>  1  aod  s  2 

2T 

s  1 

IT 

occur  anywhere  in  the  update  cycle.  These  momton  as  they  are  called,  use  the  same  obspan  as  the 
previous  run.  They  are  routinely  accomplished  when  time  permits  and  serve  to  improve  the 
prediction  process  as  a  whole.  Due  to  the  fact  that  monitor  runs  are  and  will  continue  to  be  an 
integral  part  of  the  TIP  decay  process,  the  results  of  this  research  project  will  not  be  based  solely  on 
the  mmn  TIP  decay  updates,  but  rather,  will  reQea  a  more  realistic  update  process  in  which  periodic 
monitor  runs  may  serve  as  intermediate  updates.  The  existence  of  monitor  runs  should  not  pose  a 
problem  and  is  mentioned  only  for  completeness. 

As  previously  mentioned,  one  of  the  purposes  of  TIP  processing  is  to  ensure  that  the  USSR  does 
not  mistake  a  TIP  object  for  an  incoming  re-entry  vehicle.  Consequently,  the  SSC  is  required  to 
prepare  an  OPREP-3  report  for  release  up-channel  if  a  TIP  object’s  6-hour  or  3-hour  ground  trace 
indicates  that  its  impact  point  (with  a  ±15  minute  window)  is  within  100  nautical  miles  of  the  Soviet 
border  (8,11,12,15).  Hi^er  authorities  will  release  the  information  to  the  Soviets  upon  their 
request. 
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A 


Multiple  Lmear  Recession 

Multiple  lmear  regression  is  a  method  used  to  describe  the  relationship  between  several 
mdependent  variables  and  a  dependent  variable.  It  can  also  be  used  to  predict  the  outcome  of  the 
dependent  variable  when  the  mdependent  variables  are  varied  (4:213, 18:27>31). 

If  the  mdependent  and  dependent  variables  can  be  fit  to  a  model  of  the  form  shown  in 
Equation  2,  we  can  say  that  we  have  a  linear  statistical  model  for  the  expected  value  of  the 
dependent  variable,  E(Y),  where  Y  is  the  dependent  variable.  Note  that  this,  in  turn,  means  that 
E(Y)  is  a  lmear  function  of  the  unknown  parameters  (see  Equation  2)  and  not  necessarily  a  linear 
function  of  the  mdependent  variables  themselves  (X,)  (17  495*497;.  This,  in  turn,  means  that  the 
independent  variables  could  be  used  to  predict  the  dependent  variable. 

£(  n  -  Po  *  Pi*.  *  Pj*2  *  Pj*3  ♦  PA  ♦  Ps*5  ♦  P6*6  *  PA 

where 

X,  are  the  independent  variables 

E(Y)  b  the  expected  value  of  the  dependent  variable,  Y 

P,  are  the  y-intercept  and  coeRidenta  to  be  determined 

For  thb  research  project,  the  mdependent  variables  will  correspond  to  the  early  TIP  decay 
predictions  (e,g.,  7-day  run,  4-day  run),  and  the  dependent  variable  will  correspond  to  the  final  decay 
prediction. 
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III.  Data  Coiltction  and  Descnption 


Data  Collection 

The  satellite  decay  data  used  in  this  study  was  obtained  from  the  Space  Control  Operations 
Division  (J3SOS)  at  Cheyenne  Mounts  Air  Force  Base  which  oversees  the  day-tO‘day  processing  of 
sateUite  data  by  the  Space  Surveillance  Center  (SSC).  As  previously  mentioned,  one  of  the 
responsibilities  of  the  SSC  is  to  provide  Tracking  and  Impact  Prediction  (TIP)  for  those  objects  that 
are  vvithin  fifteen  days  of  decay  and  which  have  a  high  probabUity  of  impacting  the  earth.  As  SSC 
personnel  process  each  TIP  object,  they  meticuiou^y  log  the  results  in  a  TIF  Folder  After  the 
object  htis  decayed,  the  folders  are  reviewed  by  other  SSC  crew  members  and  then  by  J3SOS 
personnel  t )  ensure  their  accuracy.  Several  kerns  from  each  TIP  folder  are  then  transferred  to 
miaofiche  for  historical  purposes.  The  items  transferred  include  the  TIP  Required  Items  Checklist, 
the  Decay  History,  the  Fmai  TIP  Message,  (he  Fmal  Element  Set,  the  Pre*/Post>£phemerides,  and 
the  Final  Ground  Trace. 

Data  Descnption 

The  data  relative  to  this  study  was  extracted  from  the  TIP  Required  Items  Checklist,  the  Decay 
History,  and  the  Hnal  Element  Set  for  each  TIP  object  studied. 

TIP  Required  Items  Checklist  The  TIP  Required  Items  Checklist  is  a  manualty'kept, 
chronolo^cal  account  of  each  TIP’s  processing.  Each  TIP  run  is  recorded  to  include  the  time  it  was 
completed,  the  aew  member  responsible,  (he  resulting  decay  prediction,  and  any  unusual 
circumstances. 

Decay  History.  The  Decay  History  is  a  computer-generated  log  of  each  special  perturbations 
update  and  includes  the  run  time,  time  of  last  observation,  epoch  time,  epoch  revolution  number,  B- 
term  value,  period,  and  resulting  decay  prediction.  An  example  is  provided  at  Figure  4. 

B-Term.  As  a  decaying  object  enters  the  atmosphere,  its  otUer  surface  interacts  with  the 
molecules  in  the  atmosphere.  This  interaction  produces  a  drag  on  the  ot^ect  uhich  can  alter  the 
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Figure  4.  Example  Decay  HUcory  (Reproduced  from  14) 


object's  speed  and  direction  as  it  descends  throu^  the  atmosphere.  The  predicted  decay  times  and 
iocationSf  therefore,  will  change  according  to  its  descent  parameters.  The  drag  on  an  object 
becomes  more  and  more  pronounced  as  it  descends  due  to  the  maease  in  thickness  of  the 
atmosphere.  The  si/e,  shape,  and  orientation  of  an  object  will  determme  its  drag  which  in  turn  wiU 
affect  us  impact  point.  As  each  object  is  updated  vnth  an  SP  MANDC,  a  term  known  as  its  B-tem 
LI  also  updated  (see  Equation  3). 


where 


B^um 


mters’ 

^biogram 


(3) 


Q  K,  the  coefficeot  of  drag 
A  is  the  relevant  surface  area 
m  is  the  mass. 

The  B'tenn  is  measured  in  meters  squared  per  kilogram  and  is  essenually  the  effective  surface 
area  per  unit  mass  affected  by  atmospheric  interaction.  The  B-ierm  is  determined  by  fittmg  C^jA/m 
and  the  >‘e!ocity  vector  to  the  observational  data.  Because  a  decaying  object’s  onentation  or 
rotatlontd  slate  can  suddenly  change,  and  because  it  may  also  br^k  apart  and  lose  pieces,  the  B> 
term  values  can  change  dramatically  throughout  the  decay  proc^ 

>Dtcay  Date  and  Time.  The  decay  date  and  time  is  legged  on  the  decay  Imtory  in  the 
Greenwich  Mean  Time  format.  It  includes  the  year,  day,  hour,  minutes,  and  seconds  (see  Figure  4). 
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V 


Latitude  and  Longitude.  The  latitudes  are  recorded  as  positive  degrees  for  north  latitudes 
and  negative  degrees  for  south  latitudes.  The  longitudes  a^c  recorded  as  degrees  east  longitude. 

Final  Element  Set.  The  Final  Element  Set  is  a  lisr'jig  of  the  parameters  n^ch  describe  the  final 
orbit  of  the  TIP  object  (see  F^ure  5).  The  data  used  in  this  study  included  both  the  final 
eccentridt)’  value  (ECC)  and  the  final  mean  motion  (XNO),  where  the  mean  motion  is 
defined  as  the  number  of  orbital  revolutions  that  the  satellite  completes  in  one  day. 


ELEKEKT  TRAMSMISSICM  FROStAM  •  OPS  ELaffiTT  DISPLAY 
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Flguru  5.  Example  Final  Elcmenr  Set  (Reproduced  from  14) 

Only  those  TIP  (Ejects  vdiich  underwent  the  entire  update  cycle  (Tniay  run  throu^  fmal  run) 
were  used.  This  was  necessary  in  order  to  accurately  analyze  the  effects  of  each  successive  update 
and  prediction.  In  the  cases  where  the  TIP  Required  Items  Checklist  and  the  Decay  History 
disagreed,  the  results  ffom  the  Decay  History  were  used.  The  data  indudes  TIP  objects  which 
decayed  during  the  years  1987  through  1990. 
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IV.  Methodology 


Final  Prediction 

The  first  objective  of  this  research  project  was  to  compare  the  accura^  of  the  TIP  decay 
predictions  made  during  1987*1990  to  the  accurate  level  asserted  by  the  Space  Surveillance  Center. 
To  achieve  this,  the  data  from  the  final  decay  predictions  was  used  as  the  control  by  which  to 
compare  the  earlier  predictions.  TUs  was  necessary  because  only  a  relatively  few  TIP  objects  are 
actually  sighted  as  they  enter  the  lower  atmosphere  or  impact  the  surface  of  the  eartL  Therefore,  a 
comprehensive  comparison  to  aaual  impact  data  was  not  possible. 

The  final  prediction  was  chosen  as  the  control  because  it  uses  observations  which  aie  closest  to 
the  actual  impact  point  and  is  considered  to  be  (he  most  accurate  prediction  available.  In  order  to 
further  justify  the  use  of  the  final  prediction  as  the  control,  a  anafysis  was  also  performed 

to  directly  compare  the  few  sighted  re-entry  pc^ts,  called  Vis  Ohs,  with  the  corresponding  final 
predictions  made  by  the  Space  Surveillance  Center.  The  Vis  Ob  data,  v^en  available,  was  taken 
from  the  Decay  History  printout  where  it  was  recorded  in  the  same  format  as  the  TIP  prediction 
data  The  specific  results  of  the  analyris  can  be  found  in  Chapter  V. 

Time  Error  Calculation 

In  order  to  assess  the  accuracy  of  the  decay  predictions,  it  was  cscesyuy  to  determine  the  amt 
error  for  each  separate  TIP  rru,  where  the  time  error  is  calculated  as  the  difference  between  the 
predicted  decay  time  for  that  TIP  run  and  the  final  run.  To  achieve  this,  it  was  first  necessary  to 
convert  the  decay  dates  and  times  (two  separate  values)  to  a  ring!#,  value.  Here,  the  hours,  minutes, 
and  seconds  were  converted  to  their  fractional  values  of  24  hours.  Toe  appropriate  day  of  the  year 
was  then  added  to  this  number  to  create  a  unique  value  representing  the  predicted  decay  date  and 
time.  These  values  were  then  used  to  calculate  a  time  error  for  eadt  run  by  subtracting  the 
predicted  decay  time  from  that  run’s  predicted  decay  time.  Note  that  predicted  decay  times  riirfifT 
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than  the  final  predicted  decay  time  yield  a  negative  time  error.  The  results  of  the  time  error 
calculations  and  corresponding  graphs  can  be  found  in  Chapter  V. 

Location  Error  Calculation 

It  was  decided  that  it  would  be  of  interest  to  also  calculate  the  difference  between  the  predicted 
location  point  for  each  run  and  the  final  prediaed  location  point  {location  error).  To  achieve  this, 
three  different  methods  were  e^qplored. 

Method  /.  The  first  method  eiqplored  attempted  to  use  the  predicted  latitude  and  lon^tude 
values  in  conjunction  wth  their  predicted  revolution  numbers  to  calculate  the  location  error.  For 
those  TIP  runs  in  vidiich  the  revoludon  number  differed  from  the  final  revolution  number,  the 
locadon  error  was  calculated  based  solely  on  the  difference  in  revoludon  number  and  the  dioaacft 
traveled  in  a  typical  decay  orbit.  Any  difference  in  locadon  within  the  same  revoludon  was 
considered  insignificant  vdien  compared  to  the  laige  location  error  incurred  by  differing  revolution 
numbeis.  It  was  then  assumed  that  the  satellite's  altitude  woiUd  not  significantly  add  to  the  size  of 
the  orbit  and,  thus,  that  the  satellite’s  orbit  could  be  considered  to  be  essentially  at  the  surface  of  the 
earth.  Therefore,  the  distance  from  the  center  of  the  earth  to  the  catellite  (r)  was  taken  to  equal  the 
mean  radius  of  the  earth  (R).  The  satellite  v^>uld  then  travel  a  distance  equal  to  the  circumference 
of  the  earth  during  each  revoludon.  A  mean  earth  radius  (R)  of  6370.949  kilometers  was  used  for  a 
calculated  earth  circumference  of  40029.853  idlometers.  (23:F'155). 

For  those  TIP  runs  in  «4iicb  the  revoludon  number  did  not  differ  from  the  fmal  revoludon 
number,  the  locadon  error  was  calculated  based  on  the  great  circle  distance  between  two  points  on 
the  earth.  Equation  4  was  used  to  calculate  the  distance  from  one  latitude  and  longitude  posidon  to 
another. 

/>  -  11142  •  [5frf(Lj)5w(Zj)  +  CosiL^)Cos{l^)Co3{X^-K{)]  {kilometers)  (4) 
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where 

D  is  the  distance 
X.  is  the  latitude 
X  is  the  longitude 

Equation  4  assumes  a  spherical  earth  and  calculates  the  number  of  degrees  along  an  arc  length 
between  the  two  points  as  measured  from  the  center  of  the  earth.  The  resulting  value  is  then 
multiplied  by  111.12  kilometers  per  degree.  Note  that  southern  latitudes  and  eastern  longitudes 
must  be  entered  as  negative  numbers  (13:72). 

This  method  was  ultimately  rejected  because  the  relative  positions  of  the  predicted  location 
points  and  the  direction  of  the  satellite’s  motion  could  not  be  accurately  determined.  Thus,  because 
a  positive  location  error  would  always  result  from  the  use  of  the  above  equation,  an  unacceptable 
bias  was  mtroduced  into  the  data. 

Method  II.  A  second  method  was  then  used  in  an  attempt  to  obtain  a  6rst  approiimation  for 
the  location  error.  Hrst,  each  TIP  otjea  was  assumed  to  be  in  a  perfectly  circular  orbit  again  at  the 
surface  of  the  earth.  Thus,  the  semi-major  ass  distance  (a)  was  taken  to  equal  the  mean  radius  of 
the  earth  (R).  Second,  each  TIP  object's  orbit  was  assumed  to  have  an  88-minute  period.  Again, 
using  a  mean  earth  radius  (R)  of  6370.949  kilometers  and  an  earth  drcumferencs  of  40029.8S3 
kilometers,  a  satellite  velocity  of  7.318  kilometers  per  second  was  calculated  (23:F-1S3).  This 
velocity  was  then  multiplied  by  the  previously  calculated  time  error  (converted  to  seconds)  to 
determine  the  distance  the  satellite  would  travel  during  that  time  period. 

Medtod  m.  The  actual  method  used  was  a  refinement  of  Method  n.  First,  in  order  to  more 
accurately  determine  the  velocity  of  each  TIP  oiqcct,  the  satellite’s  mean  motion  (n)  was  used. 
Because  data  on  the  mean  motion  for  eadt  TIP  run  was  unavailable,  an  approidmatiott  had  to  be 
made.  The  only  data  available  was  the  final  mean  motion  (XNO)  recorded  in  the  Final  Element  Set 
printout  induded  m  the  TIP  folder.  In  using  this  value  as  the  mean  motion  for  each  run,  an  error 
was  introduced  into  the  calculation  of  the  location  error.  Because  the  final  mean  motion  value  will 
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be  higher  than  earlier  values  (the  object  is  traveling  faster),  the  calculated  location  enor  will  be 
larger  than  the  actual  location  error.  In  this  vs/ay,  the  calculated  location  errors  will  serve  as  upper 
bound  error  values.  Also,  because  calculated  location  error  values  are  again  very  large  (thousands  of 
kilometers),  the  error  introduced  by  tiding  the  Hoal  mean  motion  instead  of  the  actual  mean  motion 
should  be  insignificant  vi^en  compared  to  the  location  error  values  themselves.  Any  trends, 
therefore,  that  may  exist  in  actuality  should  still  be  visible  with  the  small  additional  error  introduced 
by  strictly  using  the  final  mean  motion. 

To  calculate  the  location  error,  the  final  mean  motion  was  first  converted  to  radians  per  second 
and  then  used  to  solve  for  the  semi^major  axis  distance  (a)  u.<ing  Equation  5  and  a  value  of  3.986012 
X  10*  kmVsec^  for  n  (3:185,429). 


C 


(kilometers) 


where 

a  is  the  semi^major  ans 
^  is  the  earth's  gravitational  parameter 
n  is  the  mean  motion 


(5) 


The  satellite’s  velocity  (v)  was  then  calculated  fiom  the  semi^major  ms  Equation  6  and 
the  assumption  of  a  circular  orbit  (le.,  i  »  a)  (3:16,28).  A  circular  orbit  was  assumed  based  on  the 
mean  eccentricity  value  o!  the  TIP  objects  as  obtained  from  the  final  element  set  It  should  also  be 
noted  that  this  equation  assumes  a  two-body  motio*;  which  is  essentially  the  case  for  decaying  objects 
in  low-earth  orbit 

Each  satellite’s  velocity  was  then  multiplied  by  the  time  error  previously  calculated  (agtun 


converted  to  seconds)  to  determine  the  location  error  from  the  distance  each  TIP  object  traveled 
during  that  time  period.  Tbe  results  of  the  location  error  cJculation.41  and  the  corresponding  graphs 


v^ere 

is  the  earth’s  gravitational  parameter 

r  is  the  distance  from  the  center  of  the  earth  to  the 
satellite 

a  is  the  semi>major  axis 

^ccttftso’  Comp<trison 

To  compare  the  accuracy  of  the  1987*1990  decay  predictions  to  the  accora^  level  reported  by 
the  Space  Surveillance  Center,  it  was  dedded  to  use  the  same  accuracy  assessment  method  employed 
by  the  SSC.  This  would  allow  for  a  direct  comparison.  The  SSC  asserts  an  accura^  of  ±20  percent 
of  the  amount  of  time  until  the  object  decays.  The  “time  until  decay"  is  calculated  as  the  time 
between  that  TIP  run’s  epoch  time  (obtained  from  the  decay  history)  and  the  predicted  decay  time 
(1,9).  For  example,  if  the  amount  of  time  until  decay  for  a  particular  TIP  run  were  exactly  10  days 
and  the  time  error  calculated  was  exactly  -t-l  day,  then  the  decay  prediction  would  be  +10  percent  of 
the  time  until  decay.  The  results  of  the  analysis  for  time  error  can  be  found  in  Chapter  V. 

OPREP-3  Report  Initiation 

The  second  objective  of  this  research  projea  was  to  determine  if  it  would  be  advantageous  for 
the  SSC  to  initiate  its  OPREP-3  report  earlier  than  is  currently  done. 

The  current  procedure  for  initiation  of  the  OPREP-3  report  is  to  use  the  satellite  ground  trace 
produced  during  the  6*bour  TIP  run.  The  ground  trace  shows  the  predicted  impact  point  as  well  as 
a  ±15  minute  decay  window.  The  SSC  aew  commander  wiD  initiate  the  OPREP-3  report  if  any 
portion  of  that  ground  trace  within  the  decay  wmdow  falls  within  100  nautical  miles  of  the  Soviet 
border.  Then,  at  the  3-hour  TIP  run,  the  OPREP-3  report  is  cither  initiated,  updated,  or  cancelled 
according  to  the  updated  information  given  in  the  3-hour  ground  trace  (8,11,15). 
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To  determioe  whether  it  would  be  advantageous  to  initiate  the  report  earlier  than  at  the  6*bour 
point,  the  mean  time  error  for  the  1987*1990  time  period  was  analyzed.  The  specific  results  can  be 
found  in  Chapter  V.  Based  on  these  results,  it  was  decided  that  it  would  be  inappropriate  to  base 
the  imtiation  decision  on  data  generated  by  earlier  TIP  runs.  It  was  then  deaded  that  multiple 
linear  regression  would  be  used  to  determioe  if  a  linear  model  could  be  developed  which  would 
eliminate  biases  found  in  the  prediction  data  and  predict  the  final  decay  time  with  a  better  accuracy 
than  is  currently  obtained  by  the  6*hour  TIP  run.  It  was  surmised  that  if  a  better  accuracy  could  be 
obtained  earlier  in  the  prediction  process  using  a  linear  model  to  approximate  the  final  decay 
prediction,  the  OPREP-3  report  could  be  imtiated  earlier.  It  was  also  conjectured  that  the  existence 
of  such  a  model  might  provide  some  insight  into  the  decay  process  itself  and  suggest  an  avenue  of 
study  in  the  event  that  there  is  a  future  attempt  to  improve  the  decay  prediction  program. 

The  first  step  was  to  determine  if  a  model  of  the  form  shown  in  Equation  7  could  be  found  such 
that  the  early  decay  prediction  data  could  be  used  to  approximate  the  final  decay  prediction  time. 
The  results  are  ^ven  in  Chapter  V. 


£(</)  =  Po*  Pi't  *  P2f2*  Pj*5+  P«»4-^  f*sh* 


where 

t,  are  the  early  TIP  decay  predictions 

E(tf)  is  the  mq>ected  value  of  the  final  decay  prediction  time 

tf  is  the  final  decay  prediction  time 

arc  the  y-intercept  and  coefficients  to  be  determined 

Assuming  that  such  a  linear  regression  model  could  be  found,  the  next  step  would  be  to 
determine  which  of  the  independent  variables  could  be  eliminated  from  the  regression  model  and 
still  approximate  the  final  decay  prediction  time  with  an  accura^  similar  to  (or  better  than)  that 
obtained  currentfy  at  the  6*hour  TIP  run. 
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To  achieve  this,  a  total  of  six  linear  models,  Ei(tf)  through  EgCtf),  were  developed  to  calculate 
the  expected  value  of  the  fmal  decay  prediction  time,  E(tf),  where  each  subsequent  model 
systematically  incorporated  more  of  the  early  prediction  data  into  the  model 

The  six  models  were  then  used  to  calculate  E(tf)  for  all  180  TIP  objeas.  The  difference 
between  these  approximations  to  the  final  predktion  and  the  actual  final  predictions,  or 
approximation  error  was  then  calculated.  The  results  can  be  found  in  Chapter  V. 
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V.  Results  of  the  Analysis 


Analysis  of  Final  TIP  Decay  Predictions 

Of  the  total  180  TIP  objects  from  1987*1990  studied  in  this  research  project,  93  were  visibly 
spotted  during  reentry.  The  data  obtained  from  these  sightings  is  recorded  on  the  decay  histor>'  for 
these  objects  and  is  refened  to  as  Hr  Obs. 

The  Hnal  TIP  decay  predictions  made  by  the  Space  Surveillance  Center  (SSC)  for  thcae  93 
objects  were  compared  to  the  Vis  Ob  data.  The  mean  time  error  calculated  was  -0.87  minutes  with 
a  standard  deviation  of  12.27  minutes.  Using  Method  III  described  in  Chapter  IV,  a  location  error 
of  '301.89  kilometers  with  a  standard  deviation  of  5864.47  kilometers  was  calculated.  Note  that  there 
is  a  tendency  for  the  final  predicted  decay  time  to  be  slightly  earlier  than  the  Vis  Ob  decay  time. 

The  final  prediction  time  error  calculations  were  then  broken  down  by  year  and  are  shown  in 
Figures  6  and  7. 


MEAN  TIME  ERROR 
(Final  Run  vs  VIS  OBs) 
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FINAL  RUN  TIME  ERROR 
(Standard  Deviation) 


Flgnre  7.  Final  Time  Error  Standard  Deviation 

As  seen  in  Figure  6,  the  mean  time  error  becomes  more  and  more  negative  from  1987  to  1990. 
This  seems  to  indicate  a  growing  tendency  for  the  final  decay  prediction  time  to  be  slightly  earlier 
than  the  Vis  Ob  decay  time.  There  is  no  clear  explanation  for  this  trend.  It  may,  however,  be 
related  to  how  the  Vis  Ob  decay  point  data  is  chosen.  Figure  7  indicates  that  the  size  of  the  time 
error  deviation  deaeases  from  1987  to  1990.  A  seardi  for  a  possible  explanation  led  to  an  analysis 
of  the  level  of  solar  activity  during  the  1987'1990  time  period. 

It  was  found  that  the  level  of  solar  activity  began  to  increase  dramaticaify  in  1987  and  continued 
to  inaease  through  the  solar  maximum  which  is  believed  to  have  occurred  around  March  of  1990 
(21:3201).  Hgures  8  and  9  depict  the  sunspot  number  and  10.7  centiffleter  solar  flux  levels  during 
the  1987*1990  time  period  where  the  data  was  first  smoothed  using  a  simple  avera^ng  technique 
(22).  At  first,  this  information  would  seem  to  suggest  that  the  mean  time  error  and  standard 
deviation  should  have  increased  steadily  throughout  the  1987*1990  period.  A  closer  look  at  the  level 
of  sunspot  activity  and  solar  flux  during  that  period,  however,  revealed  that  although  there  was  an 
overall  inaease  in  both  parameters  from  1987  to  1990,  there  was  a  muck  more  dramatic  maease  in 
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those  levels  daring  1987  and  1988  than  during  1989  and  1990  (see  Figures  8  and  9).  The  trend  seen 
in  Figure  7  may  be  related  to  this  apparent  levellng-off  of  solar  activity  in  1989  and  1990. 

Again,  because  Yis  Ob  data  is  not  available  for  all  TIP  objects,  and  because  the  above  mean  and 
standard  deviation  time  error  values  were  small,  it  was  decided  that  the  fiaal  TIP  decay  prediction 
could  be  used  as  the  control  against  which  to  compare  the  earlier  TIP  decay  predictions. 

Anafysis  of  Earfy  TIP  Decay  Predictions 

Time  Errw.  The  difference  between  each  of  the  early  predicted  decay  times  and  the 
corresponding  final  prediaed  decay  time  was  then  calculated  and  plotted  for  the  1987*1990  time 
frame.  The  results  zsc  listed  in  Table  4  and  plotted  in  Figures  10  and  11. 


TABLE  4 

TIME  ERROR  RESULTS  FOR  1987*1990 


(minutes) 

wrimm 

KTUO 

636:767  1 

55r733 

2S357  1 

-jWM 

5I3I7 

mm  1 

bsiSSli 

9033 

s:767  1 

As  can  be  seen  in  Hgure  10,  there  is  a  definite  tendency  for  predicting  the  decay  time  late 
relative  to  the  final  prediction  for  all  but  the  3*hour  run.  A  possible  explanation  for  this  may  be  that 
the  atmospheric  effects  on  the  TIP  objects  are  not  being  accurately  modeled  and  that  the  overall 
drag  on  the  ot^ecc  may  actually  be  higher  than  the  current  values  calculated  by  the  SSC.  Whether 
or  not  a  sufficient  explanation  for  this  positive  bias  is  found,  however,  it  may  be  poss&ile  to 
incorporate  an  algorithm  into  the  software  used  by  the  SSC  to  model  out  the  above  bias  and,  thus, 
improve  the  overall  decay  prediction.  Figure  11  shows  a  general  decrease  in  the  time  error  deviation 
from  the  7*day  run  to  the  S'bour  im  As  expected,  this  decrease  indicates  that  the  decay  predictions 
improve  as  the  object  approadies  its  reentry  time. 
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The  prc<licted  decay  time  errors  were  then  brcVxn  down  by  year.  The  results  are  shown  in 
Figures  12  and  13.  Figure  12  again  shows  that  a  positive  error  bias  exists  for  the  decay  time 
predictions.  Figures  12  and  13  also  show  that  the  mean  time  error  for  1988  is  consistently  better 
than  for  any  other  year,  but  that  its  standard  deviation  is  consistently  higher  than  for  eny  other  year. 

A  second  look  at  Figures  8  and  9  suggest  a  possible  explanation  for  the  h^er  standard 
deviation.  As  seen  in  the  graphs,  the  level  of  solar  activity  began  to  gradually  inaease  in  1987. 
There  was,  however,  a  sharper  increase  in  those  levels  in  1988  and  then  a  leveling>off  of  those  levels 
in  1989  and  1990.  It  is  possible  that  the  larger  deviations  associated  with  the  1988  TIP  decay  dtUa 
(as  indicated  by  the  larger  standard  deviation)  may  be  associated  with  the  sharper  increase  in  solar 
activity  during  that  year  han  for  the  otner  years. 


MEAN  TIME  ERROR 
(By  Year) 


TIP  RUN 

Figure  12.  Mean  Time  Error  (By  Yeiir) 


In  an  attempt  to  explain  the  better  mean  time  error  found  for  1988,  it  was  pcstubted  that 
perhaps  a  software  fix  wa»  incorporated  by  the  SSC  to  temporarily  adjust  for  the  increasing  level  of 
solar  acthity.  However,  it  was  determined  Uu^  no  such  software  fix  was  implemented  (2).  It  was 
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TIME  ERROR 
(Standard  Deviation) 


HP  RUN 

Figure  Time  Error  Standard  Deviation  (By  Year) 


also  surmised  that  perhaps  a  procedural  change  was  implemented  by  the  SSC  trews.  This,  too,  was 
found  to  not  be  the  case  (10).  It  is  not  clear  >^y  the  mean  time  error  for  1988  is  coaisteotl^'  better 
than  for  any  other  year. 

Location  Error.  Similarly,  the  location  error  for  the  1987' 1990  time  period  was  calculated  using 
Method  III  as  descnbed  in  Chapter  IV.  The  results  are  listed  in  Table  5  and  graphed  in  Figures  14 
and  lb  Note  that  a  similar  positive  error  bias  is  seen  in  Figure  14  for  locarioa  enor  as  in  10 
for  time  error.  Again,  this  may  be  related  to  an  underestimated  drag  term  resulting  Irom  hi^  levels 
of  solar  activity.  As  with  the  time  error  calculations,  the  location  error  calculations  were  then 
broken  down  by  year.  The  results  are  shown  in  Figures  16  and  17.  As  in  Figure  13  for  the  time 
error,  it  can  be  seen  in  Figure  17  that  the  standard  deviation  for  the  locadon  error  is  also  highest  for 
the  year  1988.  Again,  this  may  be  related  to  a  sharper  increase  in  solar  activity  during  1988  than  for 
the  other  years. 
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TABLES 


LOCATION  ERROR  RESULTS  FOR  1987-1990 


1  7.Day 

1,  4-Dav 

76,'XC^ 

300,133-17 

J-E)ay 

I33,47li5~ 

l*Day 

ssisras 

12-Houf 

15325 

39^iiffi 

19^33.66 

-170 

MEAN  LOCATION  ERROR 
(1987-1990) 


Accuracy  Comparison 

As  mentioned  in  the  previous  chapter,  the  Space  Surveillance  Center  reports  a  decay  prediction 
accuracy  of  t20  percent  of  the  amount  of  time  left  until  the  TIP  object  decays.  Accordingiy,  the 
time  error  for  each  decay  prediction  was  calculated  as  a  percent  of  the  time  left  until  decay  for  that 
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LOCATION  ERROR  (1987-1990) 
(Standard  Deviation) 


MEAN  LOCAnON  ERROR 
(By  Year) 


Fi(iin  It.  Mean  Location  Error  (Ey  Year) 
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LjOCATION  error 
(Standard  Deviation) 


HP  RUN 

Figure  17.  Ix>catioQ  Error  Standard  Deviation  (By  Year) 


TIP  run.  Plots  were  created  showing  the  time  error  as  a  percent  of  the  time  until  decay  for  each 
TIP  nuL  The  results  for  the  7>day,  I'day,  and  5*hour  runs  are  depicted  in  Hgures  18'20. 

The  percentage  of  those  TIP  objects  which  actually  fell  outside  the  reported  t20  percent 
accuracy  standard  were  thee  calculated  and  plotted.  The  results  for  the  entire  1987-1990  time  period 
are  shown  in  Figure  21.  Figure  22  shows  the  results  broken  down  by  year. 

As  can  be  seen  in  Figure  22,  the  year  with  the  largest  percentage  of  its  decay  predictions  outside 
the  reported  ±20  percent  accuracy  standard  is  1988.  1990»  on  the  other  hand,  had  decay  predictions 
outside  the  reported  ±20  percent  accuracy  standard  for  only  the  2-day  and  6-bour  runs.  All  other 
decay  predictions  were  better  than  the  reported  ±20  percent  acctiracy  standard.  Again,  the  results 
for  1988  may  be  related  to  the  sharper  increase  in  solar  activity  during  th^  year  than  in  the  other 
years  (see  Hgures  8  and  9). 
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7-DAY  TIME  ERROR 


1-DAY  TIME  ERROR 
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3-HR  TIME  ERROR 
(Percent  of  Time  Until  Decay) 


Figure  20.  3-Hour  Time  Error 


DECAY  PREDICnON 
(Relative  to  20  Percent  Standard) 


DECAY  FREDICnONO  (1907-1M0) 

Figure  21.  Decay  Predioicm  Accuracy  (1987-1990) 


DECAY  PREDICTION 
(Relative  to  20  Percent  Standard) 


D 

O 


Figure  22.  Decay  Prediction  Accuracy  (By  Year) 


OPREP’3  Report  Initiation 

As  discussed  m  Chapter  IV.  the  results  of  the  accuracy  calculations  for  time  error  were  analyzed 
to  determine  bow  far  back  in  the  TIP  decay  process  one  could  go  without  significantly  lowering  the 
prediction  accuracy  used  to  initiate  the  report. 

Figure  10  suggests  that  if  the  current  TIP  prediction  process  were  employed,  it  would  be 
detrimental  to  use  the  prediction  data  generated  earlier  than  at  the  6*hour  point  to  initiate  the 
OPREP-3  report.  Use  of  the  12«hour  pre<fiction  data,  for  example,  would  incur  an  average  inacase 
in  time  error  of  12  minutes.  Such  an  increase  m  time  error  would  be  unacceptable. 

As  discussed  in  Chapter  IV,  multiple  linear  regression  was  then  used  to  first  determine  if  a 
model  of  the  form  shown  in  Equation  7  could  be  found  to  predict  the  final  decay  time.  The  results 
were  an  R-square  value  of  1 0000  and  a  p*vaiue  of  .0001.  This  means  that  at  a  significance  level  of 
05  there  erists  a  perfect  linear  relation  bet^^n  some  of  the  independent  variables  and  the 
dependent  variable  where  at  least  two  of  the  ^  tenns  are  not  zero.  The  variance  inflation  values 


37 


were  all  extremely  large,  indicating  the  independent  variables  were  all  highly  correlated  and  that  a 
gresu  deal  of  redundancy  exists  in  the  data. 

Because  it  was  found  that  a  perfect  linear  relationship  existed  between  the  early  decay 
predictions  and  the  final  decay  prediction,  six  separate  linear  models,  Hi(tf)  >  Eg(tf),  were  developed 
in  an  attempt  to  appronmate  the  fmal  decay  prediction  with  a  greater  accuracy  than  that  obtmned  b> 
the  current  TIP  decay  process  alone.  It  was  surmised  that  the  resulting  regression  model  might  be 
used  to  eliminate  some  of  the  positive  bias  found  in  the  prediction  data  and,  therefore,  improve  the 
decay  prediction  accuracy. 

Accordingly,  sbe  models  were  developed.  They  are  given  in  Equations  8  through  13.  The  first 
model  uses  only  the  7-day  prediction  data  to  calculate  the  eiqiected  value  of  the  decay 
prediction  time,  E(tf).  Each  subsequent  model  incorporates  one  additional  decay  prediction  data 
point  to  calculate  E(tf).  The  six  models  were  then  used  to  calculate  E(tf)  for  all  180  TIP  objects. 
The  difference  between  these  approximations  to  ihe  final  prediction  and  the  actual  final  predictions 
(approximaUon  error)  was  then  calculated  for  each  model. 

-  -0,U6442  ♦  0  999064(r,)  (8) 

£^{r^  -  -0.135478»0.49007(t,)*0.95U97(^)  (9) 

Eyitf  -  -0082444  +  0.41692(1, )*0.0500010i)  +  0908343(tp  (10) 

-  -0.053030 ♦  0 0Cr7471(f,) - O.Ol4395(tj) ♦  0.196518(^  .... 

♦  0,8l0492(r<)  '  ^ 

-  -0.022322  -  0.004222(f,)  -  0  005502(^  +  0.046458(t,) 

-0049000(f4)+1012311(fp  ^  ' 
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Efi)  -  -0.008370  ♦  0.000553(r,)  -  0.0029i4(^  ♦  0.0U341(^ 

~0009(m.t^ - 0.183413(c^ ♦  1.180378«<)  '  ^ 

where 

E,(tf)  is  the  expected  value  of  the  final  decay  prediction  time 

tf  is  the  final  decay  prediction  time 

tj  is  the  7<‘day  prediction  time 

t2  is  the  4-day  prediction  time 

t3  is  the  2-day  prediction  time 

t4  is  the  1-day  prediction  time 

tj  is  the  t2-hour  prediction  time 

tg  U  the  6-hour  prediction  time 

By  comparing  Figure  23  to  Figure  10,  it  can  be  seen  that  every  model  yielded  a  better  mean 
approximadon  error  than  the  TIP  runs  themselves.  For  example,  at  the  6-hour  point.  Figure  10 
shows  that  by  strictly  using  the  current  decay  prediction  procedure,  a  mean  time  error  of  9.13 
minutes  is  incurred.  Figure  23,  however,  shows  that  this  mean  time  enor  could  be  reduced  to  only 
19  seconds  if  Model  6,  E^(^),  were  incorporated  into  the  prediction  process. 

Figure  23  also  shows  that  the  best  approximation  for  the  final  decay  time  was  achieved  uring 
Model  4,  E4(tf),  where  the  estimate  of  the  mean  approximation  error  for  the  1987-1990  TIP  decay 
objects  was  zero.  Since  Model  4  uses  the  decay  predictions  generated  by  the  7-day  through  1-day 
TIP  runs,  this  would  sugge^  that  Model  4  could  possibfy  be  used  to  predict  the  &*al  decay  time  and 
decide  the  necesrity  of  OPREF-3  report  initiation  at  the  1-day  point  rather  at  the  6-hour  point 
This  in  turn  would  give  the  proper  authorities  an  additional  18  hours  notice  of  possible  satellite 
impact  near  or  within  the  Soviet  border.  Further  comparison  of  Fgure  24  with  Figure  11,  however, 
shows  that  a  greater  error  deviation  would  be  incurred  if  the  decision  were  moved  up  to  the  1-day 
point  (i.e.,  172.92  minutes  instead  of  42.23  minutes).  Assuming  an  error  Hmh  of  one  standard 


39 


deviatioa,  it  would  be  necessary  to  keep  the  decision  at  the  6’hour  point  in  order  to  main^^in 
approximately  the  same  error  limits  (i.e^  40  23  minutes  instead  of  42.23  minutes). 

Because  Figures  11  and  24  both  represent  the  amount  of  error  deviation  from  a  mean,  they 
could  be  used  to  estimate  an  error  window.  This  led  to  a  further  analysis  of  the  procedure  involved 
in  OPREP*3  report  initiation.  As  mentioned  previously,  the  current  procedure  is  to  look  at  a  ±15 
minute  window  on  the  6*hour  or  3*hour  ground  trace  to  dcteimine  the  necessity  of  the  report.  It  is 
assumed  that  this  window  was  intended  to  allow  for  some  error  in  the  predicted  decay  time 
calculations.  If  this  is  actually  the  case,  then  based  on  the  standard  deviation  results  for  the  d^hour 
point  and  using  one  standard  deviation  as  an  error  mar^n,  the  current  ±15  minute  error  window  is 
too  small.  A  more  reasonable  error  vdndow  might  be  ±42.23  minutes  for  the  current  procedure  (see 
Figure  11,  d-hour  point)  or  ±40.25  minutes  if  Model  6  is  incorporated  into  the  procedure  (see  Figure 
24). 

MEAN  APPRO;aMATlON  ERROR 
(1987-1990) 


REGRESSION  MODEL 

Figare  23.  Mean  Approximation  Error  (1987*1990) 
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VI.  Conclusions  and  Recommendations 


TIP  Decay  Accuracy 

According  to  the  results  of  the  data  used  in  this  research  project,  the  accuracy  of  the  TIP  decay 
predictions  is  pretty  much  as  reported  by  the  Space  Surveillance  Center  (Le.,  within  ±20  percent  of 
the  time  until  decay).  Except  for  the  relatively  few  predictions  which  fell  outside  of  this  margin 
(which  might  be  attributed  to  noise  within  the  data  itself),  the  decay  predictions  in  general  were 
much  better  than  the  reported  ±20  percent  as  can  be  seen  in  Figures  18-20. 

It  should  be  noted,  however,  that  a  positive  bias  seems  to  exist  mdicating  that  the  early  TIP 
decay  predictions  are  routinely  late  relative  to  the  fmal  decay  prediction.  It  may  be  possible  to 
develop  an  algorithm  to  completely  model  this  bias  out  and,  thus,  temporarily  improve  the  decay 
predictions  until  a  reasonable  explanation  and  fix  can  be  found  for  this  bias.  The  existence  of  such  a 
bias  may  be  due  to  an  increase  in  solar  activity  from  1987  to  1990  which  would  cause  an  expansion 
of  the  atmosphere  and  a  subsequent  increase  in  atmospheric  drag  on  decaying  satellites. 

Based  on  the  multiple  linear  regression  analysis  findings  of  this  research  project,  the  use  of 
linear  models  such  as  those  given  in  Chapter  V  in  conjunction  with  the  data  generated  by  the 
current  TIP  decay  process  would  allow  the  SSC  to  better  predict  the  Gnal  decay  time  by  eliminating 
some  of  the  positive  bias  found  in  the  data. 

It  may  also  be  beneficial  to  pursue  a  study  of  the  current  Special  Perturbations  model  used  by 
the  Space  Surveillance  Center  and  attempt  to  better  account  for  the  level  of  solar  activity  which 
affects  the  earth's  atmosphere  and  its  orbiting  satellites. 

OPREP-3  Report  Initiation 

With  regard  to  OPREP-3  report  initiation,  it  is  first  recommended  that  the  initiation  decision 
continue  to  be  made  at  the  6-hour  point,  rather  than  any  earlier.  It  is  also  recommended,  however, 
that  a  linear  regression  model,  such  as  Model  6,  be  used  in  conjunction  with  the  7-day  through  6- 
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hour  data  to  improve  the  Soal  decay  predictioa  time.  By  implementing  such  a  change  in  SSC 
procedures,  there  would  be  an  average  improvement  in  the  6-hour  decay  prediction  time  error  of 
5288  seconds.  Lastly,  it  is  also  recommended  that  the  tl5  minute  error  window  cunently  used  as  a 
guide  to  initiate  the  OPREP-3  report  be  widened  to  approximately  ±40  minutes  to  better  account  for 
possible  error  in  the  decay  prediction  calculations. 


43 


Biblio^phy 


I.  Barker,  BUI,  Deputy  Astrodynamics  Project  Manager.  Personal  Interview.  KAMAN  Sdences 
Corporation,  Colorado  Springs,  CO,  13  September  1991. 

1  Barker,  Bill,  Deputy  Astrodynamics  Project  Manager.  Personal  Interview.  KAMAN  Sciences 
Corporation,  Colorado  Spring.,  CO,  S  Oct(^r  1991. 

3.  Bate,  Roger  R.  cl  aL  FundaminUtls  of  Astrocfynamics.  New  York:  Dover  Publications,  Inc.,  1971. 

4  Cody,  Ronald  P.  and  Jeffrey  K.  Smith.  Appbed  Statistics  and  the  SAS  Pro^ammmg  Language 
(Third  Edition).  Amsterdam:  Elseder  ^eoce  Publishing  Co.,  Inc.,  1991. 

5  Department  of  the  Air  Force.  Operattonal  Checklist.  USSPACECOM/J3  Form  76.  TIP  Support: 
Intermediate  Runs.  OPR  J3SOS,  23  October  1989. 

6.  Department  of  the  Air  Force.  Orbitai  Anafyst  Handbook.  Colorado  Springs,  Colorado:  HQ 
AFSPACECOM,  March  1989. 

7.  Department  of  the  Air  Force.  Space  Handbook.  Maxwell,  Alabama:  Air  University,  January 
1985. 

8.  GUni,  Major  Vincent  C.,  Commander,  Space  Surveillance  Center.  Personal  interview. 

HQ  USSPACECOM/J3SO,  Cheyenne  Mountain  AFB  CO,  1  May  1991. 

9.  Hall,  Gary  L.,  Chief  Orbital  Analyst,  Span  Control  Operations  Branch.  Personal  interview. 

HQ  USSPACECOM/J3SOS,  Cheyenne  Mountain  ATO  CO,  9  September  1991. 

IQ.  Hall,  Gary  L.,  Cluef  Orbitai  Analyst,  Space  Control  Operations  Branch.  Personal  interview. 

HQ  USSPACECOM/J3SOS,  Cheyenne  Mountain  Afi  CO,  8  October  1991. 

II.  Halpin,  Lt  Col  Michael  P.,  Former  Commander,  Space  SurveUlance  Center.  Personal  interview. 
HQ  USSPACECOM/J3SS,  Peterson  AFB  CO,  1  May  1991. 

IZ  Halpin,  Lt  Col  Michael  P.,  Former  Commander,  Space  Surveillance  Center.  Personal  interview. 
HQ  USSPACECOM/J3SS,  Peterson  AFB  CO,  23  May  1991. 

13.  Hevdett-Padcard.  HP’25  Applications  iboptmis.  00025'90011,  Rev.  C  8/75,  1975. 

14.  HQ  USSPACECOM/J3SOS.  Historical  Tracking  and  Impact  Prediction  Data.  Microfiche. 
January  1987  -  December  1990. 

15.  Koch,  Lt  Col  Fred  H.,  Former  Commander,  Space  Surveillance  Center.  Personal  interview. 

HQ  USSPACEqOM/J3S£,  Peterson  AFB  CO,  I  May  1991. 

16.  Liu,  JJ.F.  "AdvaDces  in  Ortrit  Theory  for  an  Artificial  Satellite  with  Drag,”  Journal  of  the 
Asuvnautical  Sciences,  32:  165-188  (April-Junc  1983). 


44 


17.  Mendenhall,  >^^lHam  et  al.  Mathemadcal  Statistics  with  Applications  (Fourth  Edition).  Boston: 
PWS-KHNT  Publishing  Company,  1990. 

18.  Neter,  John  et  aL  Applied  Lmear  Statistical  Models:  Recession,  Analysis  of  Vanance,  and 
Ejqterimentai  Desiffis  (Third  Edition).  Boston:  Richard  D.  Irwin,  Inc.,  1990. 

19.  North  American  Aerospace  Defense  Command.  NORAD  Technical  PublicatUm,  TP  SCC  008. 
Colorado  Springs:  HQ  NORAD,  6  April  1^2. 

20.  Roehrich,  Michelle.  Orbit  PropagaU^  Types.  J^OT-TN-91-1.  Colorado  Springs: 
USSPACECOM/J3SOT,  29  August  1991. 

21.  Sawyer,  Kathy.  "Solar  Tantrum  Erupting,"  Applied  Opucs,  29:  3201-3202  (1  August  1990). 

22.  USAFETAC.  USAFETAC  Project  910m,SESS  Data  for  SaUthSe  Orbits.  USAFETAC/DNE 
Letter  vrith  two  Attachments,  15  Augu^  1991. 

23.  Weast,  Robert  C.  (editor)  CRC  Handbook  of  Chemisay  and  Physics  (Sixty-Third  Edition).  Boca 
Raton:  CRC  Press,  1982. 


45 


ytta 

Capt^  Susanne  V.  Lefebvrc  was  bom  oa  26  June  1963  in  i  Gennany.  She  graduated 

from  Fort  Walton  Beach  High  School  in  May  1%1  and  received  a  foar*year  ROTC  scholarship.  She 
attended  the  University  of  Florida,  graduating  with  high  honors  and  receiving  a  Bachelor  of  Science 
In  M^ematics  in  May  of  1996.  While  at  Che  University  of  Florida,  she  was  a  member  of  the  Air 
Force  ROTC  Detachment  150,  received  the  AFROTC  Commandant's  Award  in  May  1983,  was 
chosen  as  the  first  female  vnng  commander'in  1984,  was  a  distinguished  graduate,  and  was  awarded 
a  regular  commission  upon  graduation.  Her  fint  assignment  was  at  Lowry  AFB,  Colorado  v^re 
she  graduated  as  the  Top  Graduate  of  the  fi»t  Undergraduate  Space  Training  Class  in  February 
1987.  Upon  graduation,  she  was  then  assigned  to  the  United  States  Space  Command  from  February 
1987  to  May  1990.  During  this  period,  she  iwrked  as  an  Orbital  Analyst  in  the  Space  Surveillance 
Center,  as  the  Near-Earth  Analyst  in  J3SOS,  and  as  Deputy  Launch  OfEcer  in  J3SOS.  Her 
responsibilities  included  performing  timely  processing  of  all  tracking  data  from  a  worldwide  network 
of  26  spacetrack  sensors,  maintaining  all  procedures  concerning  near-earth  procesring,  and 
developing,  coordinating,  and  dluieminating  all  information  concerned  with  domestic  and  cooperative 
space  launches.  In  May  1990,  she  entered  the  School  of  Enguseering.  Air  Force  Institute  of 
Technology  at  Wright-Patterson  Air  Force  Base,  Ohio. 
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